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A.  M.  Geoffrion* 

A.  B.  Nelson 

Consultant  to  The  RAND  Corporation,  Santa  Monica,  California 


Numerous  requests  have  been  received  for  copies  of  the  experimental 
code  used  to  obtain  the  computational  experience  reported  in  Refs.  2 
and  3.  It  should  be  recognized  that  this  is  not  a  production  code.  It 
was  developed  to  test  the  usefulness  of  certain  innovations  applied  to 
a  simple  Balasian  algorithm.  The  central  concern  was  the  rate  of  in¬ 
crease  of  solution  time  as  a  function  of  the  number  of  variables,  rather 
than  how  to  achieve  the  smallest  possible  execution  time  for  particular 
problems.  For  this  reason,  the  simplest  possible  Balasian  algorithm  was 
used  as  the  starting  point,  and  concessions  were  freely  made  to  program¬ 
ming  expediency  (e.g.,  no  machine  language).  It  would  not  be  difficult 
to  reduce  execution  times  substantially  by  reprogramming  and  introducing 
some  of  the  more  sophisticated  tests  already  available  in  the  literature. 

We  discuss  input  in  Sec.  1;  output  in  Sec.  2;  and  give  an  example 
in  Appendix  A,  and  a  program  listing  in  Appendix  B.  For  an  outline  of 
the  working  details  of  the  algorithm,  sec  [l]  and  [2],  Familiarity  with 
these  papers  is  presumed  here. 

The  program  solves  integer  linear  programs  of  the  form 

(P)  Minimize  cx  subject  to  b  +  Ax  St  0 

Xj  ■  0  or  1 

* 

Any  views  expressed  in  this  paper  are  those  of  the  author.  They 
should  not  be  interpreted  as  reflecting  the  views  of  The  RAND  Corporation 
or  the  official  opinion  or  policy  of  any  of  its  governmental  or  private 
research  sponsors.  Papers  are  reproduced  by  The  RAND  Corporation  as  a 
courtesy  to  members  of  its  staff. 
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where  c  and  x  are  n-vectors,  b  is  an  m-vector,  and  A  is  m  by  n.  Any 
bounded  integer  linear  program  can  be  written  in  this  form,  using 
elementary  manipulations  if  necessary. 
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1 .  INPUT 


The  following  parameter  and  data  cards  appear  for  each  problem  to 
be  run: 

(a)  Parameter  card 

(b)  S-card(s) 

(c)  C-card(s) 

(d)  B-card(s) 

(e)  A-card(s) 

(f)  Blank  card. 

Problems  can  be  stacked  by  repetition  of  cards  a  through  f. 

Parameter  Card 

The  input  parameters  are: 

M  The  number  of  constraints 

N  The  number  of  variables 

L  The  number  of  variables  in  the  initial  partial 

solution  (L  must  correspond  to  the  number  of 
entries  on  the  S-card)  .  If  L  ■  0,  the  initial 
partial  solution  is  empty.  If  L  <  0,  the  initial 
partial  solution  consists  of  all  variables  fixed 
at  the  value  0. 

SC  Punch  0  if  no  imbedded  linear  program  is  desired 

(the  algorithm  then  reduces  to  a  simple  Balasian 
algorithm),  and  1  if  the  imbedded  linear  program 
is  to  be  used. 

KENUM  When  intermediate  output  is  used  (N0P  ■  0),  the 

fraction  of  all  2n  possible  solutions  that  have  been 
implicitly  enumerated  is  printed  out  every  KENUM 
times  that  backtracking  occurs.  KENUM  *=  20  is 
reasonable . 

ZBAR  If  an  upper  bound  z  on  the  optimal  value  of  the 
objective  function  of  (P)  is  known,  put  ZBAR  = 
z  -  led  +  0.0001,  where  led  is  the  least  common 
denominator  of  the  cost  coefficients  cj  (we  assume 
that  2  is  a  multiple  of  led).  Hence,  if  all  Cj 
(and  z)  are  integer,  put  ZBAR  -  z  -  .9999.  The 
effect  will  be  that  the  program  looks  only  for 
feasible  solutions  with  value  <  ZBAR.  If  no  upper 
bound  is  known,  put  ZBAR  =  0.  See  Remark  2  below. 
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ISCMAX  The  maximum  number  of  composite  constraints  that 
will  be  carried.  ISCMAX  =  4  is  -easonable. 

ISCFR  The  frequency  with  which  the  imbedded  linear  program 
is  used.  ISCFR  ■  0  means  that  it  will  never  be  used; 
ISCFR  =  j,  j  a  positive  integer,  means  that  it  will 
be  used  every  j*-*1  time.  ISCFR  =  1  has  proven  effec¬ 
tive,  but  frequently  a  value  of  8  or  so  is  even 
better. 

MAXC  If  equal  to  0,  nothing  will  happen.  If  equal  to  1, 
all  signs  on  the  C-  and  A-Cards  will  be  reversed 
automatically  when  these  cards  are  read  in.  This 
is  purely  a  convenience  for  manuscripting  and  key¬ 
punching  for  problems  with  a  preponderance  of  minus 
signs  in  C  and  A. 

MAXT  Terminates  the  calculations  after  MAXT  seconds. 

N0P  If  equal  1,  intermediate  output  will  be  suppressed; 

if  equal  0,  intermediate  output  will  appear. 

Normally  N0P  will  be  set  at  1. 

ZKBAR  Put  equal  to  led  (see  ZBAR)  minus  1.  Thus,  if  all 
Cj  are  integer,  put  ZKBAR  *  0.  The  effect  is  that 
tne  program  looks  only  for  feasible  solutions  with 
value  at  least  (ZKBAR  +  .99999)  less  than  the  best 
feasible  solution  currently  known;  this  doesn't 
exclude  any  optimal  solutions.  (A  solution  within 
A  of  the  optimum  can  be  found  if  desired  by  increasing 
the  above  value  of  ZKBAR  by  A.) 

H1,H2  Arbitrary  problem  identifiers. 


Remark  1:  The  program  is  currently  dimensioned  to  use  32,000  words  of 
core  in  such  a  way  that  the  following  limits  must  be  observed: 


M  +  ISCMAX  £  50 


N  £  90. 


Remark  2:  If  any  Cj  are  negative  (after  MAXC  has  changed  the  input 
signs,  if  it  has  value  1),  the  program  Internally  makes  a  trivial  change 
of  variables  to  make  such  c^  nonnegative:  replace  x^  by  y^  -  (1  -  x  ) 
if  Cj  <  0.  The  problem  is  solved  in  terms  of  the  new  variables,  and 
the  reverse  transformation  is  made  at  final  output  in  order  to  recover 
the  solution  to  the  original  problem.  ZBAR  must  be  set  at  a  value 


corresponding  to  the  transformed  problem  when  it  is  desired  to  use  a 
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known  upper  bound;  hence,  when  <  0  for  jcj,  put  ZBAR  »  i  -  led  +  .0001  + 

TJcl. 
jeJ1  j1 

The  fields  and  formats  of  the  parameter  card  are  as  follows: 


Parameter 


Column 


Format 


M 

1-3 

Integer 

N 

4-6 

It 

L 

7-9 

II 

SC 

10-12 

II 

KENUM 

13-17 

II 

ZBAR 

18-23 

E 

ISCMAX 

24-26 

Integer 

ISCFR 

27-29 

II 

MAXC 

30-32 

II 

MAXT 

33-37 

II 

N0P 

38-40 

If 

ZKBAR 

41-46 

E 

HI 

47-52 

Hollerith 

H? 

53-58 

II 

S-Card(s) 

The  algorithm  can  start  with  any  initial  partial  solution  (see 
[l]).  When  the  initial  partial  solution  is  desired  to  be  nonempty 
(L  >  0),  if  is  to  be  fixed  at  the  value  one  (zero)  then  "j"  ("-j") 
is  entered  on  the  S-card,  followed  by  "B"  when  an  underline  is  desired. 
The  S-card  la  divided  into  12  fields  of  5  columns  each:  1-5,  6-10, 

66-70.  Only  the  first  four  columns  of  each  field  are  to  be  used  except 
when  underlines  are  desired,  in  which  case  "B"  must  appear  in  the  fifth 
column  of  the  field. 

The  special  instruction  given  above  in  Remark  2  for  ZBAR,  when  a 
change  of  variables  is  made,  also  applies  here.  That  is,  the  sign  of 
+  )  or  +  JB  must  be  changed  when  Cj  <  0. 

C-Card  (>) 

The  values  of  the  Cj  must  be  entered  in  order  (negative  values  are 
permissible,  as  noted  above).  Each  card  has  six  fields  of  eleven  columns 
read  in  E-format.  The  fields  are  separated  by  an  unread  column  so  that 
the  values  of  the  c^  are  in  columns  1-11,  13-23,..., 


61-71. 
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B-Card(s) 

The  values  of  the  muse  be  entered  In  order.  The  format  is 
exactly  the  same  as  for  the  C-cards. 

A-Card(s) 

Only  nonzero  a^  need  be  entered,  and  they  may  be  entered  in  any 
order.  Each  value  is  identified  by  its  row  and  column.  There  are  four 
or  fewer  entries  on  each  of  the  "A"  cards.  Each  entry  has  a  seventeen 
column  field. 


Columns 

Fo rmat 

Row 

1-3 

Integer 

Column 

4-6 

Integer 

Value 

7-17 

E 

The  fields  are  separated  by  an  unread  column  so  that  the  matrix 
subscripts  and  values  of  the  a^  are  in  columns  1-17  ,  19-35  ,  37-53  ,  55-71. 
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2 .  OUTPUT 


The  preliminary,  intermediate,  and  final  outputs  are  as  follows. 

The  parameter,  "S",  "C",  and  "B"  cards  are  printed  in  that  order 
(six  values  to  a  line  for  the  "C"  and  "B"  cards).  Then  the  complete  A 
matrix  is  printed  (with  zeros),  row  by  row.  If  MAXC  =  1,  the  sign 
reversals  in  "C"  and  "A"  will  be  seen  to  have  occurred.  If  a  change  of 
variables  was  made  internally,  the  new  c,  b,  and  A  are  printed  out  (if 
no  change  of  variables  was  necessary,  the  identical  c,  b,  and  A  are 
printed  out  again  anyway) . 

If  N0P  =  0,  intermediate  output  is  produced  to  reveal  the  course 
of  the  calculations  -  each  feasible  solution  found,  each  new  composite 
constraint,  data  concerning  each  imbedded  linear  program,  and  a  summary 
of  progress  to  date  after  each  KENUM  "backtrackings."  Since  this 
information  is  likely  to  be  of  little  incremental  value  to  the  user  over 
the  final  output  information,  no  detailed  explanation  is  given  here. 

The  final  output  gives  the  problem  designation;  the  message 
"implicit  enumeration  complete"  or  "time  exceeded"  according  as 
termination  did  or  did  not  occur  within  MAXT  seconds;  the  total  exe¬ 
cution  time  in  seconds;  the  solution  (obj.  fc.  value  and  a  list  of 
which  variables  equal  1)  both  before  and  after  the  variable  change  (if 
no  variable  change  occurred,  these  solutions  are  identical);  and  6ome 
statistical  information  on  the  course  of  the  algoritiim,  such  as  the 
number  of  feasible  solutions  found,  the  number  of  times  the  imbedded 
linear  program  was  solved,  the  number  of  iterations,  and  the  time  at 
which  the  last  feasible  solution  was  found.  In  the  event  that  no 
feasible  solutions  were  found,  this  is  indicated  by  the  zeros  in  the 
solution  after  variable  change  and  the  statistic  "no.  feasible  solutions 
0."  In  the  event  that  the  time  limit  was  exceeded,  the  final  output  is 
preceded  by  a  brief  report  giving  the  proportion  of  all  2n  possible 
solutions  that  have  been  accounted  for  and  the  final  "state"  vector  [l], 
with  "B"  signifying  an  underline.  All  the  information  needed  to  restart 
the  calculations  is  available:  make  the  S-card  correspond  to  the  final 
state  vector  (set  L  accordingly),  and  put  ZBAR  equal  to  LEAST  Z  AFTER 
VARIABLE  CHANGE  -  led  +0.0001. 
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Appendlx  A 

EXAMPLE 


We  shall  illustrate  the  above  by  solving  Petersen's  fifth  example 

[4]. 

For  this  problem,  M  =  10  and  N  =  28.  We  shall  take  L  =  28,  SC  =  1 , 
KENUM  =  20,  ZRAR  =  0  (since  we  will  not  bother  to  determine  a  bound  on 
the  objective  function),  ISCMAX  -  4,  ISCFR  =  1,  MAXC  =  1  (since  we  wish 
to  avoid  keypunching  all  the  minus  signs  for  c  and  A),  MAXT  =  60, 

N0P  ■  1,  ZKBAR  ■  4  (since  the  least  common  denominator  of  the  c^  is  5), 
and  HI  »  PETE  5. 

The  S-cand  will  contain  the  numbers  1,  2,...,  28  (we  have  elected 
an  initial  partial  solution  with  all  variables  fixed  at  the  value  1). 

The  output  is  reproduced  below. 
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10  28 

28  1  20  0. 

6  1  1 

60  1  0.6000F  01  PETF  8 

P*  10 

N*  2  8 

1 

2 

* 

<1 

8 

6 

7 

8 

9 

1C 

1  1 

12 

13 

18 

16 

17 

18 

19 

20 

21 

22 

23 

28 

28 

26 

27 

-1.0000000 

E 

02 

-2.2000000 

E 

02 

-9.0000000 

E 

01 

-6.0000000 

E 

02 

-3.0000000 

E 

0? 

-6.0000000 

E 

02 

-2.0800000 

E 

02 

-1.2000000 

E 

02 

-1.6000000 

F 

02 

-8.8000000 

E 

02 

-6.0000000 

E 

02 

-1.6000000 

E 

02 

-1.0000000 

E 

02 

-1.3000000 

£ 

03 

-6.8000000 

E 

02 

-3.2000000 

E 

02 

-6.8000000 

E 

02 

-8.0000000 

t 

01 

-6.0000000 

l 

01 

-2.8800000 

E 

03 

-3.1000000 

E 

03 

-l. 1 000000 

E 

03 

-9.8000000 

E 

02 

-6.8000000 

E 

02 

-3.0000000 

E 

02 

-2.2000000 

E 

02 

-  2 . OOOOOOC 

E 

02 

-8.2000000 

E 

02 

9.3000000  E 

02 

1.2100000  E  03 

2.7700000  E 

07 

6.6200000  F 

07 

8.3700000  E  02 

8.7200000  E  02 

7.6000000  E 

02 

6.0000000  F  02 

6.7000000  E 

02 

6.9000000  E 

07 

-8.0000000 

E 

00 

-2.6000000 

E 

01 

-1.3000000 

E 

01 

-8.0000000 

E 

01 

-7.0000000 

E 

01 

-8.0000000 

l 

01 

-6.8000000 

E 

01 

-1.8000000 

E 

01 

-2.8000000 

F 

01 

-9.0000000 

F 

01 

-1.3000000 

E 

02 

-3.2000000 

t 

01 

-7.0000000 

E 

01 

-1.2000000 

E 

02 

-6.0000000 

E 

01 

-3.0000000 

E 

01 

-7.0000000 

F 

01 

-6.0000000 

E 

00 

-3.0000000 

E 

00 

-1 .BOOl  000 

E 

02 

-2 .2000000 

E 

02 

-  8. OOOOCOO 

E 

01 

-3.0000000 

E 

01 

-8.0000000 

E 

01 

-1.2000000 

E 

01 

-8.0000000 

r 

00 

-8.0000000 

F 

00 

-1.8000000 

E 

01 

-8.0000000 

E 

00 

-6.6000000 

01 

-  1 . 3000000 

E 

01 

-1.0000000 

E 

02 

-1.0000000 

E 

02 

-9.0000000 

f 

01 

-7.800O000 

E 

01 

-7.8000000 

E 

01 

-2.8000000 

E 

01 

-1.2000000 

F 

02 

-1.3000000 

E 

02 

-3.2000000 

E 

01 

-6.0000000 

f 

01 

-1.6000000 

E 

02 

-6.0000000 

E 

01 

-6.0000000 

E 

01 

-8.8000000 

F 

01 

-  1 .OOOOC JO 

C 

01 

-6.0000000 

E 

00 

-2.6000000 

E 

02 

-7.9000000 

F 

02 

-8.0000000 

E 

01 

-9.0000000 

E 

01 

-7.0000000 

E 

01 

-2. 7000000 

E 

01 

-1 .7000000 

E 

01 

-8.0000000 

E 

00 

-2.8000000 

E 

01 

-3.0000000 

E 

00 

-6.0000000 

E 

00 

-6.0000000 

£ 

00 

-2.0000000 

E 

01 

-2.0000000 

E 

01 

-3.0000000 

E 

01 

-8.0000000 

E 

00 

-3.0000000 

E 

00 

-1.2000000 

E 

01 

-1 .6000000 

E 

01 

-6.0000000 

E 

01 

-6.0000000 

E 

oc 

-3.0000000 

E 

00 

-2.0000000 

F 

01 

-8.0000000 

E 

00 

-0. 

-8.0000000 

E 

00 

-3.0000000 

E 

00 

-0. 

-2.0000000 

E 

01 

-3.0000000 

E 

01 

-6.0000000 

E 

01 

-1.0000000 

E 

01 

-0. 

-8.0000000 

E 

00 

-0. 

-0. 

-l. 0000000 

E 

01 

-8.0000000 

E 

00 

-9.0000000 

E 

00 

-6.0000000 

E 

00 

-6.0000000 

E 

01 

-3.0000000 

E 

01 

-6.0000000 

E 

01 

-1.6000000 

E 

01 

-8.0000000 

E 

oo 

-1 .8000000 

E 

01 

-7.6000000 

E 

01 

-6.0000000 

E 

01 

-1.6000000 

E 

01 

-1.1000000 

F 

01 

-3.0000000 

E 

01. 

-2.8000000 

f 

01 

-1.0000000 

E 

01 

-1.3000000 

E 

01 

-8.0000000 

F 

00 

-1.0000000 

E 

00 

-8.0000000 

E 

01 

-6.0000000 

E 

01 

-8.0000000 

E 

01 

-2.0000000 

E 

01 

-3.0000000 

E 

01 

-1.0000000 

E 

01 

-8.0000000 

E 

00 

-3.0000000 

E 

00 

-7.0000000 

E 

01 

-8.0000000 

F 

00 

-1.1000000 

E 

01 

-7.0000000 

F 

00 

-8.0000000 

E 

01 

-6.0000000 

E 

01 

-6.0000000 

E 

01 

-1.9000000 

E 

01 

-7.0000000 

E 

00 

-1.8000000 

E 

01 

-2.9000000 

E 

01 

-7.0000000 

E 

01 

-2.1000000 

F 

01 

-1.7000000 

E 

01 

-3.0000000 

F 

01 

-2.8000000 

E 

01 

-1.8000000 

E 

01 

-2.8000000 

E 

01 

-8.0000000 

l 

00 

-1.0000000 

E 

00 

-1.0000000 

F 

02 

-7.0000000 

E 

01 

-8.8000000 

E 

01 

-2.0000000 

E 

01 

-8.0000000 

E 

01 

-1.8000000 

E 

01 

-1 .8000000 

E 

01 

-6.0000000 

F 

00 

-2.0000000 

F 

01 

-8.0000000 

E 

00 

-1.1000000 

E 

01 

-7.0000000 

F 

00 

-8.8000000 

F 

01 

-6.0000000 

E 

01 

-6.0000000 

E 

01 

-2.1000000 

E 

01 

-9.0000000 

E 

00 

-1 .8000000 

E 

01 

-2.9000000 

E 

01 

-7.0000000 

E 

01 

-2.1000000 

E 

01 

-1.7000000 

F 

01 

-3.8000000 

F 

01 

-2.8000000 

1 

01 

-7. 0000000 

F 

01 

-2.8000000 

E 

01 

-8.0000000 

E 

00 

-2.0000000 

E 

00 

-1 .1000000 

E 

02 

-7.0000000 

E 

01 

-8.8000000 

F 

Cl 

-7.0000000 

E 

01 

-8.0000000 

F 

01 

-7.0000000 

E 

01 

-1.8000000 

F 

01 

-6.0000000 

E 

00 

-2.0000000 

E 

01 

-0. 

-0. 

-1 .0000000 

E 

00 

-1.0000000 

E 

01 

-6.0000000 

E 

00 

-1.0000000 

E 

01 

-0. 

-6.0000000 

E 

00 

-0. 

-6.0000000 

E 

00 

-3.2000000 

E 

01 

-3.0000000 

E 

00 

-0. 

-7.0000000 

E 

01 

-1.0000000 

E 

01 

-0. 

-0. 

-0. 

-0. 

-0. 

-3.0000000 

E 

01 

-1.0000000 

E 

01 

-0. 

-1.0000000 

E 

01 

-1.0000000 

E 

01 

-8.0000000 

E 

00 

-0. 

-1.0000000 

E 

01 

12- 


-3.0000000 

b 

00 

-4.0000000 

F 

00 

-5.0000000 

E 

00 

-6.0000000 

E 

00 

-1.2000000 

E 

01 

-1.0000000 

E 

01 

-1.2000000 

E 

01 

-1 .0000000 

F 

02 

-2.0000000 

E 

01 

-1.0000000 

C 

00 

-2.0000000 

E 

01 

-5.0000000 

F 

01 

-2.0000000 

E 

01 

-1.0000000 

E 

01 

-1.0000000 

E 

01 

-3.0000000 

E 

00 

-6.0000000 

F 

00 

-9.0000000 

E 

00 

-1.2000000 

E 

01 

-1.2000000 

E 

01 

-1.0000000 

E 

01 

-1.8000000 

E 

01 

-1.1000000 

E 

02 

-2.0000000 

E 

ot 

-2.0000000 

r 

00 

-4.0000000 

E 

01 

-6.0000000 

E 

01 

-2.5000000 

E 

01 

-1.5000000 

E 

01 

-1.0000000 

E 

01 

-3.0000000 

E 

00 

-8.0000000 

E 

00 

-9.0000000 

E 

00 

-1.6000000 

E 

01 

-1.5000000 

E 

01 

-1.0000000 

E 

01 

-1.8000000 

F 

01 

-1.2000000 

E 

02 

-2.0000000 

£ 

01 

-3.0000000 

E 

00 

-5.0000000 

E 

01 

-6.0000000 

E 

01 

-2.5000000 

E 

01 

-1.5000000 

E 

01 

-1.0000000 

E 

01 

1.0000000 

E 

02 

2.2000000 

e 

02 

9.0000000 

E 

01 

2.0500000 

E 

02 

1.2000000 

E 

02 

1.6000000 

E 

02 

1 . 0000000 

E 

02 

1.3000000 

f 

03 

6.5000000 

E 

02 

6.0000000 

E 

01 

2.5500000 

E 

03 

3.1000000 

E 

03 

3.0000000 

E 

02 

2.2000000 

E 

02 

2.0000000 

£ 

02 

-4.9700000  E 

02 

-7.7600000 

E 

02 

-4.5000000 

E 

01 

1.3000000  E 

01 

-9.6000000 

E 

01 

-2.1900000 

E 

02 

8.0000000 

F 

00 

2.4000000 

E 

01 

1.3000000 

E 

01 

4.5000000 

E 

01 

1.5000000 

E 

01 

2.8000000 

E 

01 

2.0000000 

E 

01 

1.2000000 

E 

02 

4.0000000 

E 

01 

3.0000000 

E 

00 

1.8000000 

E 

02 

2.2000000 

E 

02 

1.2000000 

E 

01 

5.0000000 

E 

00 

8.0000000 

E 

00 

8.0000000 

E 

00 

4.4000000 

E 

01 

1.3000000 

E 

01 

7.5000000 

E 

01 

2.4000000 

E 

01 

2.8000000 

E 

01 

4.0000000 

E 

01 

1.6000000 

E 

02 

4.0000000 

E 

01 

6.0000000 

E 

00 

2.4000000 

E 

02 

2.9000000 

E 

02 

2.7000000 

fc 

01 

1.7000000 

E 

01 

8.0000000 

E 

00 

3.0000000 

F 

00 

6.0000000 

E 

00 

4.0000000 

E 

00 

8.0000000 

E 

00 

3.0000000 

E 

00 

1.2000000 

E 

01 

3.0000000 

E 

00 

2.0000000 

E 

01 

5.0000000 

E 

00 

0. 

2.0000000 

E 

01 

3.0000000 

E 

01 

5.0000000 

F 

00 

0. 

0. 

5.0000000 

E 

00 

9.0000000 

E 

00 

6.0000000 

£ 

00 

1.6000000 

E 

01 

5.0000000 

E 

00 

1.8000000 

E 

01 

1.1000000 

E 

01 

3.0000000 

E 

01 

2.5000000 

E 

01 

1.0000000 

E 

00 

8.0000000 

E 

01 

6.0000000 

E 

01 

l. 0000000 

E 

01 

5.0000000 

E 

00 

3.0000000 

E 

00 

5.0000000 

F 

00 

1.1000000 

E 

01 

7.0000000 

E 

00 

1.9000000 

E 

01 

7.0000000 

E 

00 

1.8000000 

E 

01 

1.7000000 

E 

01 

3.0000000 

E 

01 

2.5000000 

E 

01 

1.0000000 

E 

00 

1.0000000 

E 

02 

7.0000000 

e 

01 

1.5000000 

E 

01 

1.5000000 

l 

01 

6.0000000 

E 

00 

5.0000000 

E 

00 

1.1000000 

E 

01 

7.0000000 

E 

00 

2.1000000 

E 

01 

9.0000000 

E 

00 

1  .8000000 

E 

01 

1.7000000 

E 

01 

3.5000000 

E 

01 

2.5000000 

E 

01 

’.0000000 

E 

00 

l. 1000000 

E 

02 

7.0000000 

E 

01 

-2.0000000 

E 

01 

-1.4000000 

E 

01 

-2.0000000 

E 

01 

-1.8000000 

E 

01 

-4.2000000 

E 

01 

-9.0000000 

E 

00 

-5.0000000 

E 

00 

-6.0000000 

E 

00 

-4.0000000 

E 

00 

-3.0000000 

E 

01 

-5.0000000 

E 

00 

-2.0000000 

E 

01 

-2.0000000 

01 

-3.0000000 

E 

01 

-2.9000000 

E 

01 

-2.0000000 

E 

01 

-3.0000000 

£ 

01 

-4.2000000 

E 

01 

-1.8000000 

F 

01 

-1.5000000 

E 

01 

-1.8000000 

E 

01 

-7.0000000 

E 

00 

-5.0000000 

E 

01 

-2.5000000 

E 

01 

-2.5000000 

r 

01 

-2.8000000 

c 

01 

-7.5000000 

E 

01 

-2.9000000 

E 

01 

-2.0000000 

E 

01 

-3.0000000 

E 

01 

-4.2000000 

E 

01 

-2.0000000 

E 

01 

-2.0000000 

F 

01 

-2.2000000 

E 

01 

-7.0000000 

E 

00 

-5.5000000 

E 

01 

-2.5000000 

E 

01 

-3.0000000 

F 

01 

-2.8000000 

E 

01 

4.0000000 

E 

02 

3.0000000 

E 

02 

4.0000000 

E 

02 

5.8000000 

E 

02 

4.0000000 

E 

02 

1.4000000 

E 

02 

3.2000000 

E 

02 

4.8000000 

E 

02 

8.0000000 

b 

01 

1.1000000 

E 

03 

9.5000000 

E 

02 

4.5000000 

E 

02 

5.2000000 

E 

02 

-1.8000000 

E 

02 

-2.5400000 

E 

02 

-2.4900000 

E 

02 

-2.5500000 

E 

02 

8.0000000 

E 

01 

7.0000000 

E 

01 

8.0000000 

E 

01 

9.0000000 

E 

01 

1. 3000000 

E 

02 

3.2000000 

E 

01 

3.0000000 

E 

01 

2.0000000 

E 

01 

6.0000000 

E 

00 

5.0000000 

E 

01 

3.0000000 

E 

01 

5.0000000 

E 

01 

1.8000000 

E 

01 

l. 0000000 

E 

02 

l. 0000000 

E 

02 

9.0000000 

E 

01 

1.2000000 

E 

02 

1.3000000 

E 

02 

3.2000000 

b 

01 

6.0000000 

E 

01 

5.5000000 

E 

01 

1.0000000 

E 

01 

8.0000000 

E 

01 

9.0000000 

E 

01 

7.0000000 

E 

01 

2.8oOOOOO 

E 

01 

2.0000000 

E 

01 

2.0000000 

E 

01 

3.0000000 

E 

01 

1.4000000 

F 

01 

4.0000000 

E 

01 

6.0000000 

E 

00 

0. 

5.0000000 

E 

00 

3.0000000 

b 

00 

4.0000000 

E 

01 

1.0000000 

E 

01 

0. 

l. 0000000 

E 

01 

4.0000000 

E 

01 

3.0000000 

E 

01 

4.0000000 

E 

01 

2.4000000 

E 

01 

6.0000000 

E 

01 

1.6000000 

E 

01 

I. 0000000 

E 

01 

1.3000000 

E 

01 

5.0000000 

b 

00 

5.0000000 

E 

01 

2.0000000 

E 

01 

3.0000000 

E 

01 

2.0000000 

E 

01 

5.0000000 

E 

01 

4.0000000 

E 

01 

4.0000000 

E 

01 

2.9000000 

E 

01 

7.0000000 

E 

01 

2.1000000 

E 

01 

1.5000000 

E 

01 

2.5000000 

E 

01 

5.0000000 

E 

00 

5.5000000 

E 

01 

2.0000000 

E 

01 

5.0000000 

l 

01 

2.0000000 

E 

01 

5.5000000 

E 

01 

4.0000000 

E 

01 

4.0000000 

E 

01 

2.9000000 

E 

01 

7.0000000 

E 

01 

2.1000000 

E 

01 

2.0000000 

E 

01 

2.5000000 

E 

01 

5.0000000 

E 

00 

5.5000000 

E 

01 

2.0000000 

E 

01 

5.0000000 

E 

01 

13 


2.0000000 

E 

01 

1 .6000000 

E 

01 

6.0000000 

E 

00 

7.0000000 

E 

01 

0. 

0. 

1.0000000 

E 

00 

l .0000000 

E 

01 

6.0000000 

E 

00 

1.0000000 

E 

01 

0. 

6.0000000 

E 

00 

0. 

6.0000000 

E 

00 

3.2000000 

E 

01 

3.0000000 

E 

00 

0. 

7.0000000 

E 

01 

1 .0000000 

E 

01 

0. 

0. 

0. 

0. 

0. 

3.0000000 

E 

01 

1.0000000 

E 

01 

0. 

1. 0000000 

1 

01 

1.0000000 

e 

01 

6.0000000 

E 

00 

0. 

1.0000000 

F 

01 

1.0000000 

E 

on 

6.0000000 

E 

00 

5.0000000 

E 

00 

2.0000000 

E 

01 

1.6000000 

F 

01 

2.0000000 

E 

01 

6.0000000 

E 

00 

1  .2000000 

E 

01 

1.0000000 

E 

01 

1 . 8000000 

E 

01 

6.2000000 

E 

01 

9.0000000 

1 

00 

1.2000000 

E 

01 

1.0000000 

E 

02 

2.0000000 

E 

01 

5.0000000 

E 

00 

6.0000000 

E 

00 

6.0000000 

E 

00 

l. 0000000 

E 

00 

2.0000000 

E 

01 

5.0000000 

E 

01 

3.0000000 

F 

01 

5 . 0000000 

E 

00 

2.0000000 

E 

01 

2.0000000 

E 

01 

1 .0000000 

E 

01 

1.0000000 

E 

01 

2.0000000 

E 

01 

3.0000000 

E 

00 

6.0000000 

E 

00 

9.0000000 

E 

00 

3.0000000 

E 

01 

2.9000000 

E 

01 

2.0000000 

E 

Cl 

1.2000000 

E 

01 

1.2000000 

E 

01 

1.0000000 

E 

01 

3.0000000 

E 

01 

6.2000000 

E 

01 

l. 8000000 

E 

01 

1.8000000 

E 

01 

1 .1000000 

E 

02 

7.0000000 

E 

01 

1.5000000 

E 

01 

1. 8000000 

E 

01 

7.0000000 

E 

00 

2.0000000 

E 

00 

6.0000000 

E 

01 

6.0000000 

E 

01 

5.0000000 

E 

01 

2.5000000 

E 

01 

2.5000000 

F 

01 

2 . SOOOOOO 

E 

01 

1.5000000 

E 

01 

1.0000000 

E 

01 

2.8000000 

E 

01 

3.0000000 

E 

00 

8.0000000 

E 

00 

9.0000000 

E 

00 

3.5000000 

E 

01 

2.9000000 

E 

01 

2. 0000000 

E 

01 

1.6000000 

E 

01 

1.6000000 

£ 

01 

1.0000000 

E 

01 

3.0000000 

E 

01 

6.2000000 

E 

01 

2.0000000 

E 

01 

1.8000000 

E 

01 

1.2000000 

E 

02 

2.0000000 

E 

01 

2.0000000 

E 

01 

2. 2000000 

E 

01 

7.0000000 

E 

00 

3.0000000 

E 

00 

5.0000000 

E 

01 

6.0000000 

E 

01 

5.5000000 

E 

01 

2.5000000 

E 

01 

3.0000000 

E 

01 

2.5000000 

E 

01 

1.5000000 

E 

01 

1 .0000000 

E 

01 

2.B000000 

E 

01 

-14 


pe  re  5 

(PL IC I T  ENUMERATION  COMPLETE  TOTAL  T IME*  4.371 

LEAST  Z  AFTER  VARIABLE  CHANGE  *  3.0950000  E  03 

0  0  0  A  5  6  7  8  0  10  II  12  13 

00000000  24  0000 

LEAST  Z  BEFORE  VARIABLE  CHANGE  >-1.2400000  E  04 

1  230000090000 

lb  17  18  19  20  21  22  23  0  25  26  27  28 

NO.  FEASIBLE  SOLUTIONS  24 
ZS  GE  ZBAR  5  TIMES 
CONSTRAINT  INFEASIBLE  15  TIMES 
AUGMENTATION  IMPOSSIBLE  2  TIMES 
AUGMENTATION  POSSIBLE  19  TIMES 
INTEGER  DUALS  0  TIMES 
LP  FATHOMED  2  TIMES 
LP  CALLED  21  TIMES 
NO.  ITERATIONS  95 

LAST  FEASIBLE  SOLUTION  AT  4.132  SECONDS 
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Appendix  B 


% IBFTC  RIP23J 

DIMENSION  A( 50,90), JF I  50,90) 

OINENSION  B(  1001. Cl 1001.BSI 1001 .SI  1001  .SB 1 1001 ,NS 1 100 1 , NF I 1 00 1 
DIMENSION  I  TEMP  I  A) , JTEMPI A) , ATENPI  A) ,  SMAX 1 100) ,  SMAXBI 100  1 .  T  I  1 00  I 
DIMENSION  CSI  100), HI  1001 
DIMENSION  XL  I  90) . Dl 90) .E I  90.90) 

DIMENSION  JHI 100), XXI 100),Y< 100)  ,PE  1 1 00 ) , K0I6 ) 

INTEGER  S.SMAX.SC.T 
COMMON  /8L S/MS l 90)  .2 BAR 
DATA  BCIB/6HB  / 

DATA  BLANK/6H  / 

100  00  110  1-1.90 
HI  1 1-0.0 
BII)«0. 0 

cin-o.o 

BSI  D-0.0 

SID-0 

SBI I l-BLANK 

NSID-0 

NFI I  )-0 

SMAX(I)-0 

SMAXBI I )- BLANK 

T  1 1 1-0 

00  110  J-l.SO 
AIJ, D-0.0 

jfij,  n-o.o 

110  CONTINUE 
II-O 
NCON-O 
NRE0-0 
N AUG-0 
NOPT-O 
N1D-0 
NAP-0 
NLPF-0 
NSIHP-0 
NFATH-0 
NENUM-0 
NTCE-0 
ITB-0 
IPOST-1 
I  INS-5 
C 

C  READ  A  NEW  SET  OF  DATA 
C  PARAMETER  CARD  FIRST 
C  »S*  CARO  THIRD 

C  'C'.'B'.'A*  MATRICES  FOLLOW  'S' 

C 

C  MINIMIZE  SUM  C(J)*XIJ) 

C  CONSTRAINTS  ARE  Bl  I  DSUM  AII,JI*XIJ)  GE  ZERO 

READ  9000.M.N.L .SC.RENUM.ZBAR.I  SC  MAX, I  SCFR.MAXC.MAXT, 

♦  N0P.ZKBAR.H1.H2 

9000  FORMAT  1 41 3,  15, E 6.0, 31 3. I  5 , 1 3,E6. 0 ,2 A6 ) 

PRINT  9993 

PRINT  900 l.M.N.L, SC, KENUM.ZBAR,  I  SC  MAX,  I  SCFR.MAXC.MAXT, 

•  N0P.ZKBAR.H1.H2 

9001  FORMAT  I A 1 3,  15, 1X.E 1 1. A, 31 3, 1 5, 1 3,E 1 1 . A , 1 X ,2A6 ) 

IF  1 1 1NS.EQ.0 )  I  INS-9999 

IF  IMAXT.EO.O)  MAXT-999999 
MAXT- 1000*MA  XT 


00000010 
00000020 
00000030 
OOOOOOAO 
00000050 
00000060 
00000070 
00000080 
00000090 
00000100 
000001 10 
00000  1  III 
00000130 
000001  AO 
00000150 
00000160 
00000170 
00000180 
00000190 
00000200 
00000210 
00000220 
00000230 
000002A0 
00000250 
00000260 
00000270 
00000280 
00000290 
00000300 
00000310 
00000320 
00000330 
000003A0 
00000350 
00000360 
00000370 
00000380 
00000390 
OOOOOAOO 
OOOOOAIO 
OOOOOA20 
00000A30 
OOOOOAAO 
OOOOOA50 
00000A60 
000O0A70 
OOOOOASO 
00000A90 
00000500 
00000510 
00000520 
00000530 
000005A0 
00000550 
00000560 
00000570 
00000580 
00000590 
00000600 


non 
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MO-M 
MI-MOM 
JSCFR- I SCFR 
ZKBAR-ZKBAR*. 99999 
PRINT  9010, M,N 

9010  FORMAT  (3H0M-, I 3, 2X, 2HN-, I  3 » 

PRINT  9992 

9991  FORMAT  I 1M  I 

9992  FORMAT  ( 1 HO ) 

9993  FORMAT  (1H1) 

Ll-L 

IF  IL.LE.O)  11-0 
READ  9100, HSIK), SBIK)  ),M1, Lit 
9100  FORMAT  (1MI4,A1II 

IF  IL.GE.OI  GO  TO  130 
L-N 

00  120  K* 1,N 
120  SOU— K 
130  CONTINUE 

READ  9200,<C(JI,J-1,N> 

9200  FORMAT  I6IEU. 0,1X3) 

IF  (MAXC.EQ.O)  GO  TO  141 
DO  140  J-l,N 

140  CIJ  1—  cm 

141  CONTINUE 

READ  9200, IBI  I  ),  I - 1, Ml 

200  READ  9400, ( ( I TEMPI K)  , J  TEMP ( K )  ,  A  TEMPI K) ) »Ml ,4) 
9400  FORMAT  ( 41 21 3, E 1 1.0,  IX I) 

END-0.0 
DO  250  Ml, 4 
KI-ITEMPIK) 

K J- JTEMPI X  I 

IF  (KI.EQ.O)  GO  TO  250 
IF  (KJ.EQ.OI  GO  TO  250 
KJF-NFIK I  )M 
NFIKII-KJF 
JF I K  I ,K JF  )«K  J 


IF  I  MAXC.NE.O)  ATEMPIK1—  ATEMPIXI 
AIK  I  ,KJ  I-ATEMPIK  ) 

END- 1.0 

250  CONTINUE 

IF  I  END. NE. 0.0)  GO  TO  200 
PRINT  9992 

PRINT  9500,  1 1  SI K ) , SB ( K  )  I,  Ml, LI 
9500  FORMAT  I  141 3X, 14, A 1 )) 

PRINT  9992 

PRINT  9600,  ICIJI, J-l.NI 
PRINT  9992 

PRINT  9600,  (  81  I  1,1-1, Ml 
PRINT  9992 
DO  251  I- 1 ,  M 

PRINT  9600, I A(  I,  J),J-1,NI 
PRINT  9991 

251  CONTINUE 
PRINT  9992 

ALL  DATA  READ  FOR  THIS  RUN 


00000610 
0000062C 
00000630 
00000640 
00000650 
00000660 
00000670 
000006B0 
00000690 
00000700 
00000710 
00000720 
00000730 
00000740 
00000730 
00000760 
00000770 
00000780 
00000790 
00000800 
00000810 
00000820 
00000830 
00000840 
00000850 
00000860 
00000870 
00000880 
00000890 
00000900 
00000910 
00000920 
00000930 
00000940 
00000930 
00000960 
00000970 
00000980 
00000990 
00001000 
00001010 
00001020 
00001030 
00001040 
00001050 
00001060 
00001070 
00001080 
00001090 
00001100 
00001110 
00001120 
00001130 
00001140 
00001150 
00001160 
00001170 
0000 1180 
00001190 
00001200 
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DO  255  J«1,N 

csm-cu) 

IF  (CIJ J.GE.O.O)  GO  TO  255 
C1J  »— CIJ  I 
DO  253  I-l.M 
B(1I*B( I)4A(  I, J) 

253  A(I,J)»-A(  I . J  I 
255  CONTINUE 

PRINT  9600, <C( J) ,J«1.NJ 
PRINT  9992 

PRINT  9600, (BII), I-l,H) 

PRINT  9992 
DO  260  I*  I , N 

PRINT  9600, ( All, J»,J«1,NI 
PRINT  9991 
260  CONTINUE 

9600  FORMAT  (61 2X, IRE  15. 8)1 

IF  I ZBAR.GT *0.0)  GO  TO  300 
ZBAR«0.0 
DO  275  J- I»N 
275  Z BAR« Z  8  AR  ♦C I  J ) 

300  ZS-0.0 

DO  325  Im  1 , M 
325  BSI I >«B<  I  ) 

DO  330  J-l.N 
330  NSIJ1-J 

IF  (L.EQ.OI  GO  TO  400 
DO  375  K-l.L 
Jl-SIX) 

K 1* IABSI J 1 1 
NSIK1I-0 

IF  CJ1.LE.0I  GO  TO  375 
ZS-ZSKIJ1) 

00  350  I-1,M 
350  BSI I l-BSI I >*A( I, Jll 
375  CONTINUE 
400  CONTINUE 

IF  (M0MSCMAX.GT.50)  1SCMAX-5O-M0 
Il-MO+ISCMAX 
DO  425  I-MI,  II 
NFI  I»-N 
DO  425  J* 1»N 
JFI I,J)-J 
425  CONTINUE 

CALL  OATIME  (0, ITO) 

IT1-IT0 
GO  TO  1910 
C 

C  INITIALIZATION  COMPLETE 
C 

1000  CONTINUE 

IF  (SC.EQ.O)  GO  TO  2400 
C  SURROGATE  CONSTRAINTS  GO  HERE 
JSCFR-JSCFRM 

IF  I  ISCFR.CT . JSCFR I  GO  TO  2400 
ML*N-L 

IF  (ML.LE.l)  GO  TO  2400 
JSCFR-0 

1050  00  1060  J«1,N 
1060  MS  I J )*0 


00001210 

00001220 

00001230 

00001240 

00001250 

00001260 

00001270 

00001280 

00001290 

00001300 

00001310 

00001320 

00001330 

00001340 

00001350 

00001360 

00001370 

00001380 

00001390 

00001400 

00001410 

00001420 

00001430 

00001440 

00001450 

00001460 

00001470 

00001480 

00001490 

00001500 

00001510 

00001520 

00001530 

00001540 

00001550 

00001560 

00001570 

00001580 

00001590 

00001600 

00001610 

00001620 

00001630 

00001640 

00001650 

00001660 

00001670 

00001680 

00001690 

00001700 

00001710 

00001720 

00001730 

00001740 

00001750 

00001760 

00001770 

00001780 

00001790 

00001800 
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NSIMP-NSIMP61 
IF(L.EO.O)  CO  TO  1076 
DO  1075  1*1, L 
J«IABS(S(  in 

1075  MS< J )  *-S< I) 

IF  (NOP.NE.O)  GO  TO  107 
PRINT  3 600 1  (  (  S(  K  I  ,  SB  (  K  I 

1076  CALL  SIMPLE  (II.  N.MO.A.C  J, XL ,0 , JH , XX , Y  ,  OB J  ,  E  ,NOP) 

IF  (NOP.NE.O!  GO  TO  1077 

PRINT  9600, OBJ. ZBAR 

1077  CONTINUE 
II*  II6IP0ST 

IF  (K0(  ll.E0.2l  GO  TO  3600 

IF  (KOm.EQ.4)  GO  TO  100 

IF  (KOI 1I.E0.6)  GO  TO  1500 

VLPS—OBJ 

IF  (VLPS.LE.  (-ZBAR ) )G0  TO  1699 
00  1350  1*1, N 

IF  (D(  I  ).NE.AINT(0(  I  D.AND.NSm.NE.O!  GO  TO  1500 
1350  CONTINUE 

00  1650  J-l.N 

IF  INS(J).EQ.O)  GO  TO  1650 

I*J 

L-L  ♦  1 

NS( J  )*0 

SB(L  )*BC1B 

IF  (01  I I.NE.O.O!  GO  TO  1600 
SIL  )—  J 
GO  TO  1650 
1600  SIL  I* J 

ZS-ZS6CI.il 
DO  1625  11-1, N 
1625  BS(I1I-BS(I1)6A(  II, J) 

1650  CONTINUE 
NID-NI061 
GO  TO  2320 
1699  KOI  l)-6 

1500  IF  (  ISCMAX.LE.O)  GO  TO  1599 
BMP 1-ZBAR 
00  1505  1*1, MO 
1505  BMP1*BMP16XL( I)*B( II 

IF  (ABS(BMPl-BIM) I .LE. 0.00051  GO  TO  1599 
IF  (M-MO.LT. ISCMAXJ  GO  TO  1520 
DO  1510  I-Ml.M 
B (  I  l-B<  I ♦  1 1 
BS( 1 1  *  B  S ( I ♦ 1 ) 

DO  1510  J-l.N 
1510  A( I , J  l*A( I61,J) 

M-M-l 

1520  B(M6l)>BMPl 

00  1550  J* 1,N 
Z JH  *XX(JI 

IF  (JH(Jl.GE.(-N)l  ZJH— ZJH 
IF  (JH(JI.GT.O)  ZJH*0. 

1550  A(M6l, J l-ZJH 
M-M6 1 

8S(MI*B(M) 

00  1575  K-1,L 
Kl-SIKI 

IF  (K1.LE.01  GO  TO  1575 


00001810 

00001B20 
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00001650 
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00001870 

00001880 

00001890 

00001900 

00001910 

00001920 

00001930 

00001960 

00001950 

00001960 

00001970 

00001980 

00001990 

00002000 

00002010 

00002020 

00002030 

00002060 

00002050 

00002060 

00002070 

00002080 

00002090 

00002100 

00002110 

00002120 

00002130 

00002160 

00002150 

00002160 

00002170 

00002180 

00002190 

00002200 

00002210 

00002220 

00002230 

00002260 

00002250 

00002260 

00002270 

00002280 

00002290 

00002300 

00002310 

00002320 

00002330 

00002360 

00002350 

00002360 

00002370 

00002380 

00002390 

00002600 
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8S(M)-BS(M)*A<M,K1) 

1575  CONTINUE 

IF  (NOP.NE.O )  GO  TO  1599 
PRINT  1598, M 

PRINT  9600, (AIM, J), J-l.N) ,B(H)  ,8  SC  M) 

1598  FORMAT  ( 22HOSURKOGATE  CONSTRA I  NTS, 2X , 14 ) 

1599  IF  (K0lll.EQ.6l  GO  TO  3400 
1900  GO  TO  2400 

1910  IJK-0 
1920  CONTINUE 

IF  (ZS.GE.ZBAR)  GO  TO  3100 
03  1950  I  1-1, MO 

1950  IF  IBS!  II  l.LT.O.O)  GO  TO  1980 
GO  TO  2320 
1980  CONTINUE 

00  2000  J-l.N 

IF  INSIJI.EQ.O)  GO  TO  2000 

IF  I ZS«C(J ) .LT.ZBAR )  GO  TO  2000 

NS( J 1-0 

L-L*l 

SBIL I-BCIB 
SIL  )— J 
2000  CONTINUE 
K INS-0 

'  ( IJK.EQ.l)  GO  TO  2220 
IIJK.EQ.2)  GO  TO  1000 
IJK-  1 

IF  (M.LT.Ml)  GO  TO  2025 
MSC-0 

I  1-MI 
I  2-M 

GO  TO  2050 
2025  MSC-l 
1 1-  1 
I2-M0 

2050  DO  2220  1-11,12 
C-BSI I ) 

DO  2100  J- 1,N 
IF  INSIJI.EQ.O)  GO  TO  2100 
IF  ( A< I,J I.GT.O.O)  Q-Q*A(I,J) 

2100  CONTINUE 

2110  IF  (Q. IT. 0.0)  GO  TO  3000 
K-NFI  I  ) 

00  2200  Kl-l.K 
Jl-JFI I , K  1 ) 

IF  I  NS  < J1 J.EQ.O)  GO  TO  2200 
2120  IF  (Q.GE.ABSIAt  I.Jim  GO  TO  2200 
NSI  J  11-0 
L-L  ♦  1 

SBIL )-BCIB 

IF  (A( I,J1).GT.0.G)  GO  TO  2150 
SIL  )--Jl 
GO  TO  2200 
2150  SIL  )-Jl 

ZS-ZS*CIJ1) 

00  2175  1 1- 1,M 
2175  BSI  Ill-BSI  1 1 1  -*•  A (  11, Jl) 

KINS-K1NSM 
2200  CONTINUE 

IF  IK INS.GE. I  INS  I  GO  TO  1920 


00002410 

00002420 

00002430 

00002440 

00002450 

00002460 

00002470 

00002480 

00002490 

00002500 

00002510 

00002520 

00002530 

00002540 

00002550 

00002560 

00002570 

00002580 

00002590 

00002600 

00002610 

00002620 

00002630 

00002640 

00002650 

00002660 

00002670 

00002680 

00002690 

00002700 

00002710 

00002720 

00002730 

00002740 

00002750 

00002760 

00002770 

00002780 

00002790 

00002800 

00002810 

00002820 

00002830 

00002840 

00002850 

00002860 

00002870 

00002880 

00002890 

00002900 

00002910 

00002920 

00002930 

00002940 

00002950 

00002960 

00002970 

00002980 
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00003000 
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2220 

CONTINUE 

00003010 

IF  (MSC.EQ.O)  GO  TO  2025 

00003020 

IF  (KINS.EQ.O)  GO  TO  1000 

00003030 

IJK*  2 

00003040 

CO  TO  1920 

00003050 

C  4A 

00003060 

2320 

CONTINUE 

00003070 

IF  (M.EQ.MO)  GO  TO  2340 

00003060 

DO  2325  l-Nl.M 

00003090 

B(  I  )*B(  I I+ZS-ZKBAR-ZBAR 

00003100 

2325 

BSI  I  )*BS(  I l+ZS-ZKBAR-ZBAR 

00003110 

2340 

ZBAR*ZS-ZKBAR 

00003120 

00  2350  J*1,N 

00003130 

2350 

SMAXIJI-SI J  > 

00003140 

GO  TO  3300 

00003150 

2400 

K  1*  0 

00003160 

DO  2500  J* 1*N 

00003170 

IF  (NS(J).EQ.O)  GO  TO  2500 

00003 1B0 

IF  ( ITB.EQ.O)  GO  TO  2430 

00003190 

IF  IZS*CI J l.GE.ZBAR)  GO  TO  2500 

00003200 

DO  2450  1*1, M 

00003210 

IF  ( A( U J ) .LE.O.OI  GO  TO  2450 

00003220 

IF  (BSm.GE.O.O)  GO  TO  2450 

00003230 

2430 

CONTINUE 

00003240 

K1-KU1 

00003250 

TIKI  I* J 

00003260 

GO  TO  2500 

00003270 

2450 

CONTINUE 

00003280 

2500 

CONTINUE 

00003290 

IF  (Kl.EQ.O)  GO  TO  3200 

00003300 

NAP*NAP»1 

00003310 

P— 1.0E10 

00003320 

DO  2575  K-l.Kl 

00003330 

J-TIK) 

00003340 

Pl-0.0 

00003350 

DO  2550  1*1,  M 

00003360 

P2*BS( I KAI  I,J) 

00003370 

IF  IP2.GE.0.0)  GO  TO  2550 

000033B0 

Pl*PUP2 

00003390 

2550 

CONTINUE 

00003400 

IF  (Pl.LE.PI  GO  TO  2575 

00003410 

P»P1 

00003420 

J1*J 

00003430 

2575 

CONTINUE 

00003440 

NS  I J  1  )*  0 

00003450 

1*1  ♦  1 

00003460 

SIL  )*Jl 

00003470 

ZS«ZS*CIJ1) 

00003480 

DO  2600  I* 1,M 

00003490 

2600 

BSI  I  )*BS( I ) ♦A (  I, Jt ) 

00003500 

H(L  )*H(  L  )  ♦  1 . 0 

00003510 

GO  TO  1910 

00003520 

3000 

NCON*NCON ♦ 1 

00003530 

C 

PRINT  3010,  I 

00003540 

3010 

FORMAT  (1H0,  I3.26HITH)  CONSTRAINT  I NFEASI BLE) 

00003550 

GO  TO  3500 

00003560 

3100 

NRED*NREOM 

00003570 

C 

PRINT  3110 

00003580 

3110 

FORMAT  ( 33HOZ  CANNOT  BE  REOUCEO  (ZS  GE  ZBARII 

00003590 

GO  TO  3500 

00003600 
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3200  NAUG«NAUG*1 
C  PAINT  3210 

3210  FORMAT  (25HONO  AUGMENTATION  POSSIBLE) 

GO  TO  3500 
3300  NOPT*NOPT ♦ 1 

CALL  DATIME  (0.IT3) 

IF  (NOP.NE.O)  GO  TO  3500 
PRINT  3310, ZS 

PRINT  3600,((S(K),SB(K)I,K*1,L) 

3310  FORMAT  (23H0  BETTER  SOLUTION  FOUNO , 5X.2HZ- , 1  PE  15 . 8 ) 

GO  TO  3500 
3400  NLPF-NLPF+l 
GO  TO  3500 

C  4B 

3500  CONTINUE 

NENUM«NENUM*1 

IF  (NENUM.LT .KENUM I  GO  TO  3530 
NENUM*0 
3505  CONTINUE 
ENUM-0.0 
00  3510  K«1,N 

3510  IF  (SB(K).EQ.BCIBI  ENUM-ENUM*. 5**K 
CALL  DATIME  (0,IT2) 

F.LT  1»  IT2-  ITO 
ELT2*  IT2-IT1 
IT  1*  IT2 

ELT  1* EL Tl/ 1000.0 
ELT2-ELT2/1000.0 
IF  (  IT2-IT0.LT. MAXT )  GO  TO  3515 
MAXT*-1 
GO  TO  3517 
351  CONTINUE 

IF  (NOP.NE.O)  GO  TO  3700 
3517  CONTINUE 

PRINT  3520, ENUM, EL T 1«ELT2,L 

3520  FORMAT  ( 1H0, F 10. 5, 38H  OF  THE  SOLUTIONS  HAVE  BEEN  ENUMERATED, 5X , 

•  15HTIME  IN  SECONOS,2X,5HTOTAL,F8.3,2X, 7HELAPSED,FB .3, 

*  5X.2HL-.I3) 

3530  CONTINUE 

IF  (MAXT.LT .0)  PRINT  3600, ( ( S(K) , SB ( K ) ) »K*l  , L ) 

3600  FORMAT  (  151 2X, 14, A  1  )  I 

IF  (MAXT.LT.O)  GO  TO  3738 

C  4B 

3700  NFATM-NFATHM 

3710  IF  (SB(L).EO. BLANK)  GO  TO  3900 
J« I ABS ( S( L ) ) 

NS( J  )*J 

IF  (S(L).LT.O)  GO  TO  3735 
ZS«ZS-C(J ) 

00  3725  I-l.M 
3725  BS (  I  )*BS( I )- At  I, J) 

3735  S8(L  )*BLANK 
S(L  )»0 
L«L-1 

IF  (L.GT.O)  GO  TO  3710 
C  FINISHED 

3738  CONTINUE 

PRINT  3739,H1,H2 

3739  FORMAT  (1H1,5X,2A6) 

DO  3740  J-l,N 


00003610 

00003620 

00003630 

00003640 

00003650 

00003660 

00003670 

00003680 

00003690 

00003700 

00003710 

00003720 

00003730 

00003740 

00003750 

00003760 

00003770 

00003780 

00003790 

00003800 

00003810 

00003820 

00003830 

00003840 

00003850 

00003860 

00003870 

000038B0 

00003890 

00003900 

00003910 

00003920 

00003930 

00003940 

00003950 

00003960 

00003970 

00003980 

00003990 

00004000 

00004010 

00004020 

00004030 

00004040 

00004050 

00004060 

00004070 

00004080 
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00004110 
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00004160 
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00004190 

00004200 


-23- 


3740  S(J  1*0 

00004210 

00  3742  J*1,N 

00004220 

K-IABSI SMAXI J ) ) 

00004230 

IF  CK.EQ.OI  GO  TO  3744 

00004240 

3742  SIX  1*1 

00004250 

3744  00  3  7  46  K-l.N 

00004260 

IF  ISIKl.NE.O)  GO  TO  3746 

00004270 

SMAXI J)*-K 

00004280 

J*J*1 

00004290 

3746  CONTINUE 

00004300 

CALL  OATIME  (0, IT2) 

00004310 

ELT  1*  IT2-  ITO 

00004320 

ELT  I*  EL  T 1  / 1000.0 

00004330 

IF  (MAXT.LT.O)  GO  TO  3752 

00004340 

PRINT  3750*  EL  T I 

00004350 

3750  FORMAT  ( 30H0IMPL IC I T  ENUMERATION  C OMPLE TE ,5X , 1 IHTOT AL  TIME*, F0. 31 

00004360 

GO  TO  3758 

00004370 

3752  PRINT  3755, ELTI 

00004380 

3755  FORMAT  (14H0TIME  E XCEEOEO , 5 X , 1 IHTOTAL  TIME*,F8.3) 

00004390 

3758  CONTINUE 

00004400 

ZBAR*  ZBAR *ZK BAR 

00004410 

PRINT  3760*  ZBAR 

00004420 

3760  FORMAT  I 32H0LEAST  Z  AFTER  VARIABLE  CHANGE  *,1PE15.8> 

00004430 

1*0 

00004440 

3800  DO  3810  K»i,N 

00004450 

3810  TIK  1*0 

00004460 

DO  3820  K*1,N 

00004470 

K  1*  I ABS  (  SNA XI K  )  ) 

00004480 

3820  IF  (  SMAXI K  )  . GT . 0 )  TIK1)*K1 

00004490 

PRINT  3830, ITIKI,X-1,NI 

00004500 

3830  FORMAT  I15I4X, 131) 

00004510 

IF  1  I .NE.O)  GO  TO  3845 

00004520 

Z8AR-0.0 

00004530 

DO  3835  J-l.N 

00004540 

X* I ABS I SMAXI  J  1  ) 

00004550 

IF  (CS(K).LT.O.O)  SHAXIJ)*- SMAXI J) 

00004560 

IF  ISMAXI JJ.GT.O)  ZBAR-ZBAR+CSI K) 

00004570 

3835  CONTINUE 

00004580 

PRINT  3840, ZBAR 

00004590 

3840  FORMAT  I 33HOLEAST  Z  BEFORE  VARIABLE  CHANGE  -,1PE15.8) 

00004600 

1*1 

00004610 

GO  TO  3800 

00004620 

3845  CONTINUE 

00004630 

ELT 3*  IT 3-  ITO 

00004640 

ELT3-ELT3/1000.0 

00004650 

NITER-NFATH+NFATH-l 

00004660 

PRINT  3850,N0PT,NRE0,NC0N,NAUG,NAP,NID, 

NLPF.NSIMP, NITER, ELT3 

00004670 

3850  FORMAT  (23H0N0.  FEASIBLE  SOLUTI ONS,  15/ 

00004680 

•  1 1 H  ZS  GE  ZBAR,  15, 6H  TIMES/ 

00004690 

•  22H  CONSTRAINT  INFE A SIBLE , 1 5 ,6H 

TIMES/ 

00004700 

♦  24H  AUGMENTATION  1 MPOSSI BLE , 1 5 , 

6H  TIMES/ 

00004710 

*  2?  AUGMENTATION  POSSIBLE, 15, 6H 

TIMES/ 

00004720 

•  14  INTEGER  0UALS.I5.6H  TIMES/ 

00004730 

*  12.  P  F  A  THOME  0, 1 5,6H  TIMES/ 

00004740 

•  10H  P  CALLED, 15, 6H  TIMES/ 

00004750 

*  15H  NO.  ITERATIONS, 15/ 

00004760 

*  26H  LAST  FEASIBLE  SOLUTION  AT,F8.3,9H  SECONDS  I 

00004770 

GO  TO  100 

00004780 

3900  SB  I L l-BCIB 

00004790 

SIL  )*-S(L  ) 

00004800 
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J-IABSISILI) 

IF  <  S ( L  I.GT.O)  GO  TO  3950 
ZS*ZS-C(J  ) 

DO  3925  I" l»M 
3925  BS( I)-6S( I)-A(I.J) 

GO  TO  1910 
3950  ZS-ZS+CU) 

DO  3975  I«i,M 
3975  BS! 1)-BS( IMA!  I  •  J  I 
GO  TO  1910 
END 

% IBFTC  SIMPLE 

C  AUTOMATIC  SIMPLEX  REDUNDANT  EQUATIONS  CAUSE  INFEASIBILITY 

SUBROUTINE  SIMPLE!  INFLAG, MX, NN, A,B»C ,KO,KB,P, JH,X,Y,OBJ,E,NOP) 
REAL  B(1),C!  1),P(  1),  X(l),Ytl)  .OBJ 
REAL  E( 90, 90 ) 

INTEGER  INFLAG,MX,NN,K0<6) ,KB< 1),JHI l) 

FQU I VAl r  CE  (XX, LL) 

C  THE  FOLLOWING  DIMENSION  SHOULD  BE  THE  SAME  HERE  AS  IT  IS  IN  CALLER. 
REAL  A! 50,90) 

REAL  AA,AIJT,B8«C0ST»DT,RC0ST,TEXP,TPIV»TY»X0LD»XX,XY,YI»YMAX,EH 

INTEGER  I, I A, INYC, IK, I  TER, J, JT,K,KBJ,LL,M,N, JT2 

INTEGER  NCUT ,  NUNVR ,NVER , NUMPV 

LOGICAL  TR IG, VER 

LOGICAL  F INV.FFRZ , SCH 

COMMON  /BLS/MSI90I.ZBAR 

DIMENSION  NF ( 90) 

C  SET  INITIAL  VALUES,  SET  CONSTANT  VALUES 

FINV  -  .FALSE. 

TRIG  -  .FALSE. 

ITER  -  0 
LPSEQ  »  LPSEQ+1 
NUNVR  >  0 
NUMPV  -  0 
M  »  MX 
N  -  NN 

TEXP  ■  .5**16 

NVER  ■  M/2  ♦  5 

NCUT  «  4*M  ♦  10 

IF  ( INFLAG. EQ.O)  GO  TO  1410 
C  IMPOSE  CORRECT  TEMPERATURE  ON  ROWS 

FFRZ  * .TRUE. 

00  1960  L»l,M 

IF  (MS(L).EQ.NF(L) )  GO  TO  1955 

IF(MS(L)*NF(L).GT.O.OR.(MS(L).EQ.O. AND. X(L).GE.O. ) )  GOTO  1950 
I«L 

IF  (NF(L).NE.O)  GO  TO  1925 
1920  IF  (JH(I).GT.O)  GO  TO  1930 
C  IF  JH  DISAGREES  WITH  MS  DO  SPECIAL  PIVOT 

IF  (MS(L  I.GT.O. AND. JH(L) .GE.(-M) )  GO  TO  1950 
IF  (MS! L  ) . L T .0 .AND. JH! L).LT.(-M))  GO  TO  1950 
C  SPECIAL  PIVOT,  SWITCH  SINGLETONS 


1925 

00  1926 

J"  1  »N 

P(J1  » 

PIJ)  ♦ 

El  I, Jl 

E  (  I ,  J  ) 

•  -EM, 

J) 

1926 

CONTINUE 

OBJ  -  OBJ  ♦  X!  I ) 

X(I)  «  - 

X!  1 ) 

JHL  -  JH(L  1 

IF  UHL. 

GE.(-N)) 

JH!  L ) 

00004810 

00004820 

00004830 

00004840 

00004850 

00004B60 

00004870 

00004880 

00004890 

00004900 

00004910 

00004920 

00004930 

00004940 

00004950 

00004960 

00004970 

00004980 

00004990 

00005000 

00005010 

00005020 

00005030 

00005040 

00005050 

00005060 

00005070 

00005080 

00005090 

00005100 

00005110 

00005120 

00005130 

00005140 

00005150 

00005160 

00005170 

00005180 

00005190 

00005200 

00005210 

00005220 

00005230 

00005240 

00005250 

00005260 

00005270 

00005280 

00005290 

00005300 

00005310 

00005320 

00005330 

00005340 

00005350 

00005360 

00005370 

00005380 

00005390 

00005400 
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IF  (JHL.LT.I-MH  JHILI  «  -L 

00005410 

GO  TO  1950 

00005420 

c 

00  FULL  PIVOT  ON  SINGLETON 

00005430 

1930 

JT  *  -I 

00005440 

COST  -  P(I) 

00005450 

IF  (MSI  H.GT.Ot  GO  TO  1931 

00005460 

JT  ■  JT-M 

00005470 

COST  *  1  .-COST 

00005480 

1931 

EN  -  1. 

00005490 

GO  TO  630 

00005500 

C 

GET  COLUMNUT) 

00005510 

1932 

SCH  -  .FALSE. 

00005520 

IF  ( COST »GT .0*1  GO  TO  1938 

00005530 

1935 

GO  TO  1000 

00005540 

C 

SELECT  ROW! IR) 

00005550 

1936 

IF  ( IR.NE.O.OR.SCHI  GO  TO  1940 

00005560 

SCH  ■ .TRUE . 

00005570 

1938 

EN  — EN 

00005580 

00  1937  J*1,N 

00005590 

YIJI  •  -Y( J  I 

00005600 

1937 

CONTINUE 

00005610 

GO  TO  1935 

00005620 

1940 

I F  (  (  SCH. AND. AB  S( COST I.GT.TPIVl.OR.IR.EO.O)  GO  TO  1980 

00005630 

1941 

IF  (EN.GT.O. 1  GO  TO  1945 

00005640 

00  1942  J  *  1 »M 

00005650 

YIJI  •  -YIJI 

00005660 

1942 

CONTINUE 

00005670 

1945 

GO  TO  901 

00005660 

C 

PIVOT!  1R.JT) 

00005690 

1950 

NFILI  »  MS(L I 

00005700 

1955 

IFIJHIL l.LT.OI  GO  TO  1960 

00005710 

IA-JHIL  I 

00005720 

KBI  IAI-L 

00005730 

1960 

CONTINUE 

00005740 

FFRZ  •  .FALSE. 

00005750 

GO  TO  910 

00005760 

C* 

START  WITH  SINGLETON  BASIS 

00005770 

1410 

00  1402  J-l.N 

00005780 

KBIJI  -  0 

00005790 

1402 

CONTINUE 

00005800 

FFRZ  ■. FALSE. 

00005810 

1400 

00  1401  I  -  l.M 

00005820 

JMII)  »-I 

00005830 

NFIII  -  MSI  1 1 

00005840 

IF  INF!  1 1  .LT.O.OR.  ( NF I  I )  .  E  Q.  0.  AND.  B  ( I  1 .  LT.O.  1 1  JHIII—  I-N 

00005850 

1401 

CONTINUE 

00005860 

C* 

CREATE  INVERSE  FROM  'KB •  ANO  'JH»  (STEP  71 

00005870 

1320 

VER  -  .TRUE. 

00005880 

INVC  ■  0 

00005890 

NUMVR  -  NUMVR  ♦ 1 

00005900 

TRIG  -  .FALSE. 

OOOC59 10 

OBJ  »  0. 

00005920 

00  1113  I  -  l.M 

00005930 

00  1151  J-l.N 

00005940 

El J.  1 1  -  0. 

00005950 

1151 

CONTINUE 

00005960 

IF  ( JH(  II.LT.I-MII  GO  TO  1111 

00005970 

IF  (JH( II.GT.O)  JHIII  -  0 

00005980 

Eli, 11  -  1. 

00005990 

P(  II  ■  0. 

00006000 
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x< n  *  bid 

GO  TO  1113 

nil  eii.ii  *  -l. 
pm  *  ♦  l . 

OBJ  *  OBJ  ♦  8(1) 

xm  =  -biii 

1113  CONTINUF 

00  1102  JT  =  l.N 

IF  <KB(  JT) .60.01  GO  TO  110? 

GO  TO  600 

C  GET  COLUMN ( J  T  J 

1 1 1 A  TV  r  TP  I V 
I R  *  0 

COST  *  C(JT) 

00  1104  I  =  1 » M 

COST  =  COST  ♦  A( JT, I ) *P I  I ) 

IFIJHI  I) .NE.O.OR.Xf  I ).NF.O..OR.ABS(Y(  I ) ).LF.TY)  GO  10  1 104 
TY  =  ABSI Y( I  )  ) 

IR  =  I 

1104  CONTINUE 

IF  (IR.NE.0I  GO  TO  1119 
TY  *  0. 

DO  1106  I  *  1  •  M 

IFfJHf I ) .NE.O.OR.Xf I ).EQ.0..0R.ABS(Y( I ) J.LE.TPIV)  GO  TO  1106 
IF  (ABSfYI  m.LE.TY«ABS(Xf  I  )))  GO  TO  1106 
TY  *  ABSI Yf ll/XIIII 
IR  *  I 

1105  CONTINUE 

1119  IF  ( IR.NE.O)  GO  TO  900 
C  PlVOTflR.JT) 

FINV  *  .TRUE. 

IF  (NOP.EO.O)  PRINT  1199.LPSE0 

1199  FORMAT! 15H0INVERT  FAIL  IP,  14) 

GO  TO  1410 

1102  CONTINUE 

C*  PERFORM  A  SIMPLEX  ITERATION 

1200  VER  »  .FALSE. 

500  00  503  I  •  I'M 

IF  <NF(  D.EO.O.ANO.Xf  D.LT.O.  )  XIIUO. 

503  CONTINUE 

C*  FIND  MINIMUM  REDUCED  COST  ISTEP  3) 

599  JT  x  o 

BB  =  0.0 

00  701  J  ■  1 ,  N 

IF  (KB(J).NE.O)  GO  TO  701 
OT  -  cm 
DO  303  I  «  1 , M 

OT  *  OT  ♦  A ( J, I ) *P <  I  ) 

303  CONTINUE 

IF  (DT.GE.BB)  GO  TO  701 
BB  x  OT 
JT  ■  J 
701  CONTINUE 

DO  702  I* 1 ,M 

IF  I JHf I I.LT.O)  GO  TO  702 
IF  fP(I) .L T .BB )  GO  TO  703 
IF  ( (  l.-PI  m.GE.BB)  GO  TO  702 
BB  »  l.-PI  I) 

JT  x  - | -M 
GO  TO  702 


00006010 
00006020 
00006030 
00006040 
00006050 
000060t>0 
00006070 
00006080 
00006090 
00006100 
00006110 
00006120 
00006130 
00006140 
00006150 
00006160 
00006170 
00006 1  BO 
00006190 
00006200 
00006210 
00006220 
00006230 
00006240 
00006250 
00006260 
00006270 
000062B0 
00006290 
00006300 
00006310 
00006320 
00006330 
00006340 
00006350 
00006360 
00006370 
00006380 
00006390 
00006400 
00006410 
00006420 
00006430 
00006440 
00006450 
00006460 
00006470 
00006480 
00006490 
00006500 
000065  10 
00006520 
00006530 
00006540 
00006550 
00006560 
00006570 
00006580 
00006590 
00006600 
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703 

B8*PI 1 1 

000066  i  o 

JT  *  -I 

00006620 

702 

CONT 1NUE 

00006630 

COST  =  BB 

00006660 

IF  (JT.E0.01  GO  TO  2 03 

00006650 

IF  (  ITER.GE.NCUT I  GO  TO  160 

00006660 

ITER  *  ITER  ♦  1 

00006670 

C* 

MULTIPLY  INVERSE  TIMES  A(.,JT>  (STEP 

6) 

000066BO 

IF  (JT.LT.O)  GO  TO  630 

00006690 

C 

BEGIN  SUBROUTINE  GET  COLUMN(JI) 

00006700 

600 

DO  610  1*  1  » M 

00006710 

Y I I  I  =  0.0 

00006720 

610 

CONT INUE 

00006730 

DO  605  I*  l.M 

00006760 

A  I  J  T  =  A! JT, I ) 

00006750 

IF  (AIJT.EQ.O.)  GO  TO  605 

00006760 

DO  606  J  *  l.M 

00006770 

Y ( J  )  *  Y ( J  )  ♦  A  I J  T  *E ( J, II 

000067B0 

606 

CONTINUE 

00006790 

605 

CONTINUE 

00006B00 

GO  TO  660 

00006810 

630 

JT2  *  -JT 

00006R20 

EM  *  1 . 

00006830 

IF  (JT2.LE.M)  GO  TO  631 

00006860 

JT2  *  JT2  -  M 

00006850 

EM  *  -l. 

00006860 

631 

DO  632  1*1, M 

00006870 

YIII*  EM*E(  I, JT21 

00006880 

632 

CONTINUE 

00006890 

660 

YMAX  *  0. 

00006900 

DO  620  I  *  l.M 

00006910 

YMAX  *  AMAXK  ABSI  Y(  I  II.YMAX  I 

00006920 

620 

CONTINUE 

00006930 

TPI V  *  YMAX  *  TEXP 

00006960 

C 

END  OF  GET  COLUMN 

00006950 

IF  (FFRZ)  GO  TO  1932 

00006960 

IF  ( VFR I  GO  TO  1116 

00006970 

RCOST  *  YMAX/BB 

00006980 

IF  (TRIG. AND. BB. GE.  I-TPIV) 1  GO  TO  203 

00006990 

TRIG-BB.GE.I-TPIVI 

00007000 

C* 

SELECT  PIVOT  ROM  (STEP 

5) 

00007010 

1000 

AA  *  TP  IV 

00007020 

IR  *  0 

00007030 

1002 

DO  1003  1  *  l.M 

00007060 

IF  (XII  I .NE.O..OR.YI II.LE.AA.OR.NF(I).NE.O)  GO 

TO  1003 

00007050 

AA  *  Y( I » 

00007060 

IR  =>  I 

00007070 

1003 

CONTINUE 

00007080 

IF  I IR.NE.O)  GO  TO  1020 

00007090 

AA  *  0. 

00007100 

00  1010  1  «  l.M 

00007110 

IF  INFII) .NE.O.OR.YI 1 I.LE.TPI V. OR. YIII . IE . AA*X I  I 11  GOTO  1010 

00007120 

AA  -  Y (  1  1  / XI  1  I 

00007130 

IR  =  I 

00007160 

1010 

CONT INUE 

00007150 

1020 

IF  ( FFR 2  I  GO  TO  1936 

00007160 

IF  I IR.EQ.O)  GO  TO  207 

00007170 

C* 

PIVOT  ON  (  IR, JT)  (STEP 

6) 

00007180 

901 

IA  •  JH( IR 1 

00007190 

IF  (IA.GT.O)  KB  f I A )  =  0 

00007200 

C  BEGIN  SUBROUTINE  PIVOT! IR.JT) 

900  NUMPV  *  NUMPV  ♦  1 

JH<  IR)  ■  JT 

IF  I JT .GT .0)  KB! JT  I  =  IR 
Y I  *  -Y! IR  I 
YIIR)  *  -1.0 
DO  904  J  *  1,M 
XY  *  E<  IR, J  > / Y I 
IF  (XY.EQ.O.  I  GO  TO  904 
P  IJI  •  P  ( J )  ♦  COST  *  XY 
E (  IR. J)  *  0. 

00  906  I  >  1 ,  M 

El  I,J  )  «  El  I. J)  ♦  XY  ♦  YIII 

906  CONTINUE 
904  CONTINUE 

XY  *  XI  IR  »  ✓  Y I 
DO  908  I  -  1 «  N 
XOLD  *  XII) 

XI n  *  XOLD  ♦  XY  *  Yl I  I 

908  CONTINUE 
YIIR)  *  -Yl 
XI IR )  *  -XY 

C  END  OF  PIVOT 

OBJ  *  OBJ  ♦  X Y*CO ST 
IF  (VER)  GO  TO  1102 

C  EXCHANGE  ROUS  IF  SLACK  PIVOTED  IN  WRONG  ROW 

IF  I JT.GT.0.0R.JT2.EQ.IR)  GO  TO  907 
XY  -  XI  IR  I 
XIIR)  «  XIJT2) 

XIJT2I  *  XY 
00  909  I  *  1 » M 
XY  *  El  IR,  I) 

El  IR,  I  )  *  El JT2,  I  ) 

El JT2, I )  -  XY 

909  CONTINUE 

IA  »  JHIJT2) 

JHIJT2)  -  JT 
J  HI  IR)  =  IA 
KBIIA)  *  IR 

907  1NVC  *  INVC  *1 

C  TO  STEP  1  IF  NOT  INVERTING,  TO  STEP  7  IF  INVERTING 
IF  IFFRZ)  GO  TO  1950 
IF  I  OBJ .GE  .  Z  BAR  I  GO  TO  180 
IF  I F  IN V )  GO  TO  1200 

910  IF  I  INVC.GE.NVER  )  GO  TO  1320 

GO  TO  1200 

C*  END  OF  ALGORITHM,  SET  EXIT  VALUES 

207  IF  (RCOST.LE. (-1000. > )  GO  TO  203 

C  INFINITE  SOLUTION 

K  =  2 

GO  TO  250 
180  K«6 

GO  TO  250 

C  PROBLEM  IS  CYCLING  PERHAPS 

160  K  «  4 

PRINT  161.LPSEQ 

161  FORMAT! 31  HO ITERATION  LIMIT  EXCEEDED  ON  LP.I4) 

GO  TO  250 

C  FEASIBLE  OR  INFEASIBLE  SOLUTION 

20  3  K  ■  0 


00007210 
00007220 
00007230 
00007240 
00007250 
00007260 
00007270 
00007280 
00007290 
00007300 
00007310 
00007320 
00007330 
00007340 
00007350 
00007360 
00007370 
00007380 
00007390 
00007400 
00007410 
00007420 
00007430 
00007440 
00007450 
00007460 
000074  1 0 
00007430 
00007490 
00007500 
00007510 
00007520 
00007530 
OCOO  7540 
00007550 
00007560 
00007570 
00007580 
00007590 
00007600 
00007610 
00007620 
00007630 
00007640 
00007660 
00007660 
00007670 
00007680 
00007690 
00007700 
00007710 
00007720 
00007730 
00007740 
00007760 
00007760 
00007770 
00007780 
00007790 
00007800 
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250  00  1399  J  -  1,N 
XX  «  0.0 

KBJ  -  KB( J ) 

IF  IKBJ.NE.O)  XX  -  X( KB J  I 
KB ( J  )  -  LL 

1399  CONTINUE 
KOdl  >  K 
K0( 2 )  *  ITER 
KOI  3  I  -  INVC 
KOI  4 )  »  NUMVR 
KOI  5 1  -  NUNPV 
KOI  6 1  «  JT 
IF  (NOP.NE.O)  RETURN 
PRINT  162, LPSEO.IKOU), 1-1,6) 

162  FORMAT  I 3H  LP,I5,6H  KO  ,616) 

C  PRINT  1982 

1982  FORMAT  1 21H0I  JH  NF  MS  ,P,Y,X,B/1X) 

C  00  1963  I-l.M 

C  PRINT  1984, I, JHII),NF(l>,MStI>, Pill, YU), XII), Bill 

1983  CONTINUE 

1986  FORMAT  I  IX, 413,  4F 12.6) 

RETURN 

1980  IF  INOP. EQ. 01  PRINT  1981, IPSEQ.l,! R.SCH, COST 


00007810 

00007820 

00007830 

00007840 

00007850 

00007860 

00007870 

00007880 

00007690 

00007900 

00007910 

00007920 

00007930 

00007940 

00007950 

00007960 

00007970 

00007980 

00007990 

00008000 

00008010 

00008020 

00008030 


1981  FORMAT!  3H0CP, 14, 12H  FAIL,  SLACK, 13, 4M  IR-I3.5H  SCH-L1.3H  C-F19. 6 >00008040 
IF  IIR.NE.OI  GO  TO  1941  00008050 


GO  TO  1410 


00006060 


END 


00006070 
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